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TECHNICAL MEMORANDUM 1346 


. NONLINEAR THEORY OF A HOT-WIRE ANEMOMETER™ 


By R. Betchov 


We study here the properties of a hot-wire anemometer under the 
supposition that the heat transfer from the wire to the air depends, 
first, on the difference in temperature and, second, on the square of 
that difference. This latter hypothesis is confirmed by experience, 
and the consequences might be of great importance, that effect of non- 
linearity is stronger than the effect of thermal conduction. 


I. THE NONLINEAR LAW OF KING 


The heat quantity © removed per second by an air stream V from 
& wire of the diameter d and unit length is given by King in the form 


Q = (a + b/Vd)T (1) 


where T denotes the temperature difference between wire and air. King's 
calculation, approximately confirmed by experience, yields 


a=! b = V2xx!B'c'! (2) 


with x' = thermal conductivity of the air, 6' = density of the air, 
and c! = specific heat of the air for constant volume. 


These quantities may vary with Т, and experience shows that a 
increases while b remains practically constant. Intuitively, one may 
interpret this effect by saying that the air in contact with the wire 
is heated which increases its conductivity. In compensation, its density 
decreases because the pressure varies only very slightly. Obviously, 
the effects on x' and б! compensate one another, and only a varies. 


*"Thedrie non-linéaire de l'anémométre à fil chaud." Koninklijke 
Nederlandsche Akademie van Wetenschappen. Mededeling No. 61 uit het 
Laboratorium voor Aero- en Hydrodynamica der Technische Hogeschool te 
Delft. Reprinted from Proceedings Vol. LII, No. 3, 1949, pp. 195-207. 
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King states in his original report (ref. 1) that a increases by 
0.111 percent per degree; he also describes there an effect of the 
diameter on that term a which. we shall not discuss here. Thus it is 
advisable to write 


Q = fall ут) + в/чарт (3) 


where y takes the nonlinearity into account. 


One must not forget the hypotheses on which King bases his calcu- 
lation: he contends that the air flow is without viscosity and that 
the heat flow in the immediate proximity of the wire is constant. He 
uses the specific heat at constant volume although the pressure is 
certainly more constant than the density. For that reason, we consider 
equation (3) as an empirical relation, valid for the wire unit length, 
and would wish to see King's problem made the subject of a more thorough 
investigation. 


Here we intend to study the effect of the term у on the properties 
of the hot wire; we simplify the notation by introducing P so that 


P = Руа (4) 
Q = а(1 + P+ yT)T (5) 
II. GENERAL EQUATION OF THE HOT WIRE 


We shall use the following symbols: 


S resistance of the wire, per unit length, at operating 
temperature 

So resistance of the wire, per unit length, at ambient temperature 

I intensity of the electric current heating the wire 

a coefficient of the variation of S according to the 
temperature 


m weight of the wire, per unit length 
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c specific heat of the wire, in joule/gram and degree 


хо. . coefficient of the thermal conductivity of the wire in watt/ cm 
and degree | р 


С wire cross-sectional агеа 
1 semilength of the wire 
t time 
x coordinate of position, varying from 1 to -1 
We put 
A= 58-01 + P) = nme (6) 


The equation of the hot wire must express the equilibrium between 
the heat supplied per second, the heat removed by the air stream, the 
heat required to raise the temperature of the wire, and the heat trans- 
mitted by conduction. One obtains 


ате (еа y lae me c Om LUN. DE 
a2S 2 
о 


For the steady-state case, and introducing the parameters 


2 
x ж NO А ~ І 
yo 1* = Z = ———(S - So) 
1 aso(A - 12) 5012 


(8) 








к-н. МИШ 
3 3 


а25 (1 - I2/A)2 ` 
one obtains the equation (T) in the form 


ESCH 


2+ З с22 - 2% 
ду 


1 (9) 
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III. EXACT INTEGRATION OF THE STATIC CASE 


Multiplying equation (9) by д7/ду and integrating, one obtains, 
with a constant 


92 = Vaz3 + 72 - 27, + const (10) 


ду 


Опе notes that 97 /ду is zero for а negative value of Z and may 
be zero for two positive values of Z. At the ends of the wire, one 
has 5 з 8, and 7 = О; at the center, Z must be positive and 97,/ду 
zero. The range of interest for us lies, therefore, between 7 - 0 and 


the first positive root which gives 22 = 0 which we shall denote by 
y 


Z-B. We put 


Z(y) = B - X*(y) (11) 
By virtue of the relation 
GBJ + B? - 2B + const = О (12) 
one obtains 
(ax/dy)® = -GX* + EX? + р (13) 
with 
Е = 1 + 3GB D = 2(1 - В) - 3GB2 (14) 


Following, we shall consider B as a new integration constant 
indicating the temperature in the middle of the wire. At the center of 
the wire, one has X = О and 


h(Ox/dy)° = D (15) 


"M 
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which shows that D is positive and generally small. 


the wire, Z = 0 and X = +В. 


The roots of equation (13) are 


Introducing the parameter В so that 


sinh В = sh в = 20 


one can write equation (13) in the form 


k(ax/dy)® = spen В - 1) - 5 E B-1)+ vi 


We define the angle Фф, function of у, so that 


X= D sin 中 
T VE ch Vi: е 
pee l- k^sin^ 


with 


k2 = ch B +1 
2 ch В 


Equation (18) then becomes 


At the ends of 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 
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and we obtain the elliptic integral of the first kind 


Ф 
— —À ar an ру (22) 


The variable у varies from -& to 8, with 


yee 
Ш 


1/1% (23) 


and we have, for у = Ё and х2 = В 


Í к ch = (2h) 


with 
| sin Ф 
Ув = 27: € A (25) 
Bch В 1 - квір 
This last equation may be written 
ete ea (26) 
А с 
sin“? ах 
From equation (20) one may deduce 
= = 1 + tanh2(B/2) (27) 
k 


With B ranging from O to зо, k^ varies from 1 to 0.5. 
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With б and B known, one may calculate successively D, E, В, 
k, and Gase, A table of U(k,p) then permits us to calculate Ё. 


- Figure 1 gives the results obtained by this procedure and allows - 
starting out from a prescribed wire and with y known - to determine 
B from С and E, 


The temperature distribution over the wire is given by Z asa 
function of у, or by X as а function of ф and ф а function of у. 
The relation Ф(у) is given by the quotient of equations (22) and (24), 
namely 


T ар 
y do Vi - ksin% (28) 
: m i 

o Vi - к25іп2Ф 


From equations (19) and (25), one obtains a relation between Ф and X, 
namely 


2 
X? sing l-k sin Phas 


-Sig M ex (29) 
B Sin ^9, ak 1 - к25іп2ф 
The total resistance R of the wire is given by 
+1 2 а 
2891 
r= f в ах = 2801 дж Z dy + 2852 (30) 
-1 A-I О 


We introduce the cold resistance К and the function X 


о 


RI? E 
і В - Во = ——— | -F x2 з (31) 
o 
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We replace X according to equation (19) and dy according to equa- 
tion (21), namely 


2 中 
Rt 2D JE max sin 


R - Ro = —— lB - 2 
GE E(E ch B)3/2 Jo (1 - k2sin2p)3/* E a 


The integral is equal to 20U/dk© and we obtain 


R-R 


(33) 


RoI? : ор W аю 
ни E 


A - I2 E ch В U 


The values of 3U/Sk2 can be deduced from а good table of U(k,@) 
with sufficient approximation. 


If the wire were perfect, the expression in parentheses in equa- 


tion (33) would have to be replaced by unity; therefore, we shall intro- 
duce the quantity M so that 


2 
м = 1 - в + —2D_ 90/98 (3h) 
E ch В U 


The formula (33) then gives us 





2 
l - MI^/A 
R-R (35) 
° 1- I/A 
and the important relation 
P GNE. lc ep sore AS (36) 
R - Ro l -M 


This last equation permits easy determination of the wire charac- 
teristics because R, R,, and I can be measured accurately and because 


the curves obtained as functions of R/R。， for instance, indicate directly 
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the effects of conduction and of nonlinearity. We calculated the values 


of 


г-и as а function of G and 8; figure 2 shows our results. 





A good approximation is given by 


M v d (31) 








where B, corresponds to B, in the case $= om 


Bo зУї + 6G -1 (38) 


3G 


IV. A FEW USEFUL APPROXIMATIONS 


2 


In performing the calculations necessary for the plotting of fig- 
ure 1, we have noted that one may assign to k the value unity without 
introducing large errors. 


This implies В = О and equation (19) then gives 


Х = ур tan Ф (39) 


The integral (22) becomes 


中 
U = 2 _ .l,lcsen?.-wWgy, (40) 
О 2 1 - sin 中 2 


whence, one deduces 


ch(2U) = ch(VEy) (41) 
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x? = 元 (ch VEy - 1) (42) 
At the limits, one has 
В = (ch VEE - 1) (43) 
which gives 
Z= — — l1- ch Бу (44) 
l - l/ch VEE ch VEE 


In order to calculate М, one must put 





P . : 
lim du/ae = Sch sin’? ap = 8100 - mie (45) 
k—1 2 Јо cos 4 cos“ 
which gives 
кеше Ёш d Th VEE (46) 
1 - 1/ch VEE VEE 


One can see that, due to the nonlinearity, Ё is replaced by VEE, 
and the central temperature is lowered. 


If one takes equation (44) as solution of equation (9), one sees 
that the equation is satisfied for a term of approximately (ch VEy/ch VEE)®, 
and that B is given approximately by 


В = (v1 6-1 - l/ch VEE) (47) 
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with 
RT E > V1 + бс (48) 


We shall take an example that represents an extreme case. We 
choose a platinum wire with 10 percent of iridium, a diameter of 了 microns, 
and a length of 22 = 1.14 mm. Exposed to an air stream of 5 m/sec and 
heated with 75 mA, it gives us 


Р = 2.9 А = 1.2 х 10-2  1?/A = 0.47 


1% = 0.15 mm t 


3.75 


Assuming у = 1.2 x 1073 and with the aid of figures 1 and 2, we 
determined 


G = 0.25 B = 0.76 E = 1.56 M = 0.4 


The other parameters have the following calculated values: 


k 79.50 


0.99756 Фах 


sh B = 0.10 р 


0.0177 U = 2.34 


In figure 3, we show the profile of the temperatures calculated 
exactly, the profile calculated with the approximation k = 1, and the 
profile calculated with у = 0, It can be seen that the nonlinearity 
offers a more uniform temperature distribution, and that the approxi- 
mation is sufficient. 


By means of equation (35) one calculates = = 1.5, whereas the 
о 


calculation with у = O would give the result 1.9. Тһе mean tempera- 
ture giving i- = 1.5 would be 380°. 
о 


As to the term RI^/R - Ro, it changes from the value 1.24 A when 
the current is very weak to the value 1.35 A when I attains 75 mA. 
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It varies therefore by about 10 percent between = = 1 and & = 1.5 
о 


which is of the same order as the variations observed. 


Thus this magnitude is constant at 10 percent and King's law is 
verified; however, we shall see later on that the thermal inertia is 
very different from the expected value. 


V. THE DYNAMIC EQUATION OF THE HOT WIRE 


In order to calculate the variation of the total resistance of the 
wire when the current or the air stream fluctuate, one must go back to 
the equation (7). Replacing in this equation I, 8, and V (contained 
in the term A) by І + jet, 8 + веб, апа ү + ved®t, one obtains 
after suppressing the terms of the order zero, two, and more as well as 
the factor ejot 





: 2 aP v 2ay mc хс дв 
2811 + Ifs = —— -(S - So) + As + (S - So)s + — jos - 5— 一 一 
2aS V а28,2 або або ax? 
(49) 
We introduce 
_ т? 
7 з 5 А - ТС (50) 
2 
Sar 


and, identical to the formula (9) of our Mededeling No. 55 (ref. 5) 


we = Pog I8) (51) 
mc 


One then obtains, after introduction of 2 


; 2 2 
гічної CU pil Р 2 абі + 302 + Aner) _ 9-® 
I 1 - I/A 1 - І2/А ду? 
(52) 
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The function Z(y) appears twice in this differential equation, and we 
must utilize the approximation k= 1 in order to avoid great compli- 
cations. Writing the expression Z according to equation (44), one 
obtains 


Ss 3GB M зов | m 
зу + efa + T 1/ch Vite + Anden - ch VEL - 1 ch vir) = Су + Co ch VEE 

(53) 
with | 


М 2 
набір оба іа уе = 


1-19/Al -l/chVEy 2 V|L Я Pj _ 12/, 1 - 1/ch VEE 


(54) 
Gah ИРС _ Чё. eg ПЕВНО НО Ык. SAE 
а Ty - mini - 1/са УЕРО 2 V + Р ү _ т2/д 1 - 1/ch VE 
(55) 


The solution without the second member of equation (53) is а 
Mathieu function taken for a purely imaginary value of the argument, 
but if МЕЁ is large enough, for instance, more than 4, one may neglect 
the term with ch VEy of the member on the left side of equation (53). 
Actually it is important only when у is close to Ё; however, we shall 
take as a limiting condition z(£) = О, and the product z ch VEy will 
never be important. We shall also neglect the term 1/ch МЕР in the 
denominator of one of the terms on the left side of equation (53) which 
amounts to taking Z = B in the factor term of z. 


Taking the definition of E into account, one then has 


.9'z , й ch VEy 
ay? + z(E + jo/o*) = C1 + Co ch VEE (56) 


the solution of which - null for у = ЖЕ - is 


C с 
0 б _ ch ply) e E VEy ch SC (57) 


ES Ep? ch pVEE jo/w*\ch VEE ^ ch pVEE 
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with 
p= Vl + до/Еауж (58) 


It can be seen that due to the nonlinearity the characteristic 
frequency w* is replaced by а new frequency which is higher by a fac- 
tor E. (Compare with equation (8) of ref. 5.) The amplitudes according 
to Cj and Co also are modified by the nonlinearity. 


The variation r of the total resistance R will be 


+ RoI 4 ft 
r= f в ак = i | z dy (59) 
-1 A - I? 5 ug 


The integration gives 


RoI Cy tanh pyEé Co 1 tanh vet) 
ae Oe I SS [Wu ee PVE. |+ VER. mM pyss 
r a ch E | 十 5 Al anh МЕ! (60) 


We assumed above that МЕЁ is sufficiently large; also, we may 
assign to the hyperbolic tangents the value unity (the presence of a 
complex argument is here not of importance). 


One then obtains 


з 区 (1 - 1/pVEt) -2 E (61) 
а ee — - + 一 -一 一 一 ~ 一 一 
“人 Ep? mes до/аж я р ) 


Vi. RESPONSE TO A FLUCTUATION OF THE CURRENT 


If one assumes v = О in equations (54) and (55), one may trans- 
form equation (61), suppressing the terms containing 1/ch VEE 


г _ _21ВоТ/А | (1 - Gi - в) - i/pVEt) SCI ажур - Ї 
portu is 75 
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This formula gives the alternating electromotive force produced 
at the boundaries of the hot wire by the modulation current i in 
addition to the normal electromotive force Ri. Although this formula 
seems to be complicated, it can be adapted to the needs of practice. 
If the frequency w/ 2x tends toward infinity, equation (62) becomes 


ep 21R, 1^/A (1 - I2/A(1 - В) - ІС/АВ/УКЕ (63) 


no. e/a)? An 


If one takes into account that equation (35) gives R, that equation (46) 
gives M, and that equation (51) gives w*, one finds 


rl Tb. aon? with (pue (64) (65) 
0/?n лис 


The constant С corresponds to that of our former publications, 
and equation (6h) shows that the electromotive force гї taken at high 
frequency permits to measure C without being impeded by either con- 
duction or nonlinearity. The method described in reference lh is there- 
fore indicated rather than the one consisting of measuring the phase 
displacements, with w being of the order of w*. 


One may immediately verify this point by assuming w in the equa- 
tion (19) as very large, thus making the effect of the terms of conduction 
and of nonlinearity negligible. 


When w tends toward zero, the electromotive force becomes 


т) (66) 


Thus the complex function rI according to equation (62) will change 
from equation (64) into equation (66) when w varies from O to a large 
value. The complex trace of this function gives practically a semicircle 
which permits to put approximately 


rt = EHO = О) Це = Ze (67) 
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where w** denotes the effective characteristic frequency. In fig- 

ure 4, we plotted the semicircle and indicated a few values of w/w**, 
From equation (62), and in the case of the preceding example, we calcu- 
lated the electromotive forces rI for a few values of w/Ew*. One 

can see that the two functions blend, at low frequencies, if one assumes 
ож = 1.1Е0% = 1.720%. At high frequency, equations (64) and (67) will 
be equal which permits calculation of a satisfactory value of w** when 
the effect of thermal inertia is important. One obtains 


dee = (68) 
e NETTES зе z B) 
1 - 1/VEE 
1 - І2/АМ 


In the case treated one finds w** = 1.23Eo*, that is, w** = 1.920%: 
the effective characteristic frequency is almost twice the expected 
value; therefore, the approximation (67) gives & correct plot of the 
function, but the phases according to equation (68) will be only within 
a lO-percent accuracy. 


The denominator of equation (68) depends chiefly on £, and it 
increases the characteristic frequency. Instead of compensating each 
other as in the static case, the two effects reduce the thermal inertia. 


Intuitively, one may say that the conduction shortens the hot part 
of the wire and thus reduces the heat required for modifying the central 
temperature; the nonlinearity depends on the presence of hot air around 
the wire, and the thermal inertia of the air is negligible which improves 
the spherical response of the anemometer. 


When the wire is dusty, the quantity of immobile air is greater, 
and the experience shows that the term a in equation (3) is increased 
while b remains unchanged. One must therefore expect a dynamic action 
of the dust of the wire to the extent that E is modified. The dynamic 
effect may be more important than the static effect. 


The wire in the quoted example demonstrates that, with a term 
RI*/R - Ro constant at 10 percent, the characteristic frequency may be 


almost twice the normally foreseen value. 
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VII. RESPONSE TO A FLUCTUATION OF THE AIR STREAM 


Assuming i = О in the formulas (54) and (55), and maintaining у, 
one may transform equation (61) into 


Бор ЕЕ oe 
гҮ 1 + Р (1 - I2/A)? 1 - 1/ch VEE 
(69) 
[eE cr | 
With œ tending toward zero, one has 
1-31. 
BI ^? Vs 
pj BW des А а mc B „з. — З Vee (70) 
2vVi+P (1. 19/4)? 1 - 1/ch VEE E 
If w tends toward infinity, one has 
d 
3 doe 
INC MIENNE ES __ VEE (71) 
?V1l-*P(i.712/A)? 1 - 1/ch VEE | do/o** 
and one may approximately replace equation (69) by the semicircle 
rI = гІ(о = 0) (72) 
1 + jo/ov*** 
with 
week = Sue VER. (73) 


1-21/VRE 


which, in the case of the example treated above, gives шк = 1.81o*. 
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Thus there is, on principle, no equality between the dynamic reac- 
tion to a variation і and toa variation v. 


This difficulty arises due to the term VEE, namely to the conduction; 
the nonlinearity tends toward diminishing its importance (factor VE). 


We hope to publish some empirical results, and the calculation of 
the differences indicated by different authors, in the near future. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 1.- Relation between the magnitudes В, С, and Є. 
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Figure 2.- Values of M/1 -M according to С and t. 
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Figure 4.- Study ofthe tension rl. The semicircle corresponds to 


formula (67); the points give гі, according to (62), for several values 
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